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ABSTRACT 

We present results from Chandra X-ray Observatory and Rossi X-ray Timing Explorer (RXTE) 
observations of the magnetar IE 1547.0—5408 (SGR J1550— 5418) following the source's outbursts in 
2008 October and 2009 January. During the time span of the Chandra observations, which covers days 
4 through 23 and days 2 through 16 after the 2008 and 2009 events, respectively, the source spectral 
shape remained stable in the Chandra band, while the pulsar's spin-down rate in the same span in 
2008 increased by a factor of 2.2 as measured by RXTE. This suggests decoupling between the source's 
spin-down and radiative changes, hence between the spin-down-inferred magnetic field strength and 
that inferred spectrally. The lack of spectral variation during flux decay is surprising for models of 
magnetar outbursts. We also found a strong anti-correlation between the phase-averaged flux and 
the pulsed fraction in the 2008 and 2009 Chandra data, but not in the pre-2008 measurements. We 
discuss these results in the context of the magnetar model. 

Subject headings: pulsars: individual (IE 1547.0-5408, PSR J1550-5418, SGR J1550-5418) — stars: 
neutron — X-rays: bursts 



1. INTRODUCTION 

Anomalous X-ray pulsars (AXPs) and soft gamma re- 
peaters (SGRs), though previously thought to be differ- 
ent classes of objects, are now believed to all be strongly 
magnetized neutron stars, known as "magnetars". They 
are characterized by long spin periods of 2-12 s and large 
spin-down rates that imply ultra-strong su rface dipole 
mag netic fields of 10 14 -10 15 G (see reviews bv lKaspill2007l 
and [Mercghctti 2008]). The most remarkable feature of 
magnetars is their violent outbursts, during which the X- 
ray luminosity could increase by a few orders of magni- 
tude. In the context of the twisted magnetosphere model 
( Tho mpson et a l. 2002), the energy release is due to mag- 
netic field re-arrangement, which is possibly triggered by 
crustal deformation. While the post-outburst behavior 
of a magnetar can provide important information on the 
physical conditions of the magnetosphere and the stellar 
surface, only a handful of follow-up studies have previ- 
ousl y been carried out with focusing X-ray instruments 
(e ^lWoods~et"aLl 12001; IGavriil et all [2001 llsrael et al.l 
l2007f ). because the transient nature of these events re- 
quires prompt observations. In this work, we study re- 
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cent outbursts from the AXP IE 1547.0-54080 in 2008 
and 2009 using observations made with the Chandra X- 
ray Observatory and the Rossi X-ray Timing Explorer 
{RXTE). 

The X-ray source IE 1 547.0—5408 was discovered 
by lLamb fc Ma rkert (1981) with the Einstein Observa- 
tory. Based on the X-ray spectrum and infrared flux, 
IGelfand k, Gaenslerl (|2007l ) first suggested the source as 
a magn etar candidate. Th e detection of radio pulsa- 
tions bv lCamilo et al.l (|2007|) directly confirmed the pul- 
sar nature of this object; its spin period of 2.1s is 
shorter than any other known ma gnetarsB The pul - 
sar's spin-down rate as reported bv lCamilo et al.l (2007) 
implies a surface ma gnetic dipole field of 2.2 x 10 14 G. 
lHafpern et al.l (|2008l ) reported a high state of this mag- 
netar in 2007 based on XMM-Newton observations, and 
they concluded that an X-ray outburst had occurred be- 
tween 2006 and 2007. On 2008 October 3 (MJD 54742), 
IE 1547.0—5408 s howed bursting ac tivity with outbursts 
detected by Swift (jlsrael et al.ll2010T) and by the Gamma- 
ray Burst Monitor ( GBM) onboard the Fermi Gamma- 
ray Space Telescope (|Kaneko et al.ll2010l ). On 2009 Jan- 
uary 22 (MJD 54853), the AXP entered a second active 

2 Also known as SGR J1550-5418 or PSR J1550-5418. 

3 See the online magnetar catalog at 

http : //www. physics .mcgill . ca/~ pulsar /magnetar /main, html 
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TABLE 1 

Chandra Observation Parameters and Results 



Obs. 


Date 


MJD 


ObsID 


Days Since Exposure 


Count Rate a 


Pulsed Fraction a 










Outburst 


(ks) 






2008 October 


1 


2008 Oct 7 


54746.6 


8811 


4.2 


12.1 


1.35(1) 


0.21(1) 


2 


2008 Oct 10 


54749.5 


8812 


7.2 


15.1 


1.17(1) 


0.22(1) 


3 


2008 Oct 18 


54757.6 


8813 


15.2 


10.1 


1.15(1) 


0.33(1) 


4 


2008 Oct 21 


54760.8 


10792 


18.4 


10.1 


1.07(1) 


0.35(1) 


5 


2008 Oct 26 


54765.1 


8814 


22.8 


23.0 


0.99(1) 


0.31(1) 


2009 January 


6 


2009 Jan 23 


54855.0 


10185 


2.0 


10.1 


0.95(l) b 


0.09(1) 


7 


2009 Jan 25 


54856.7 


10186 


3.7 


12.1 


3.1(2) 


0.09(1) 


8 


2009 Jan 29 


54860.8 


10187 


7.8 


13.1 


2.5(1) 


0.13(1) 


9 


2009 Feb 06 


54868.6 


10188 


15.6 


14.3 


2.2(1) 


0.12(1) 



a In the 1 — 7keV energy range. 

k Obs. 6 was made with the HETG, which results in a reduced count rate. Correcting for the effect 
gives a count rate of ~ 3.6. 
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phase. More than 200 bursts were detected within a few 
hours by Swift, International Gamma-Ray Astrophysics 
Laborat ory (INTEGRAL), F e rmi GBM, and Suzaku 
WAM (iGronwall et all 120091: iMereghetti et~aT1 120091: 
Savchcnko ct al. 2010; Kancko ct al. 2010; Terada ct al. 



20091 ). Follow-up imaging observations with Swift and 
XMM-Newton revealed dust scatte ring- X-ray rin g s cen - 
tered on the source, from which iTiengo et all (|2010f ) 
deduced a source distance of 3.9 kpc. Based on 
Suzaku obse r vation s taken 7 days after the 2009 outburst, 
lEnoto et al.l (|2010f) reported a hard power-law tail in the 
spectrum, of photon index T between 1.33 and 1.55, and 
extending up to at least 110 keV. Hard X-ray pulsations 
were also detected by INTEGRAL in the 20-1 50 keV 
band (|Kuiper et al.l [20091: Iden Hartog et al.l 12009ft . The 
pulsed emission has T — 1.55 and the spectral shape re- 
mained stable over the observation period covering from 
day 2 to day 9 after the outburst. 

In this paper, we report on results from Chandra 
and RXTE observations taken after the 2008 and 2009 
events. The observations and results are reported in Sec- 
tion [21 and we discuss their physical implications in Sec- 
tion [3] We summarize our findings in Section [4j 

2. OBSERVATIONS AND RESULTS 

2.1. RXTE Observations and Results 

Since the 2008 October outburst, IE 1547.0-5408 has 
been monitored regularly with RXTE. Data are col- 
lected with the Proportional Counter Array (PCA) in- 
strument that consists of five collimated xenon/methane 
multi-anode Proportional Counter Units (PCUs). Only 
GoodXenonwithPropane mode data were used in our 
analysis, which give 1 /is time resolution and 256 energy 
channels in the 2-60 keV energy range. We considered 
events from only the top Xenon layer of each PCU to 
maximize the signal-to- noise ratio (S/N). Typical expo- 
sure times were in the range 2-7 ks. In all, we report 
here on a total of 55 observations taken between MJDs 
54743 and 54911. A more complete report on the RXTE 
data will be presented elsewhere (R. Dib et al., in prepa- 
ration). 

2.1.1. RXTE Timing 

To determine the pulsar spin parameters, we cleaned 
the RXTE data and selected events between 2 and 



6.5 keV. As the source exhibited many short X-ray 
bursts, we removed all burst intervals for this analysis. 
The photon arrival times were first corrected to the so- 
lar system barycenter, then binned with 31.25 ms time 
resolution. The time series were folded at the nominal 
spin period and pulse arrival times were extracted by 
cross-correlating with a template. These arrival times 
were then fitted to a simple phase-coherent timing model 
using the TEMPO software packageQ Details on the 
phase fol ding and periodic ity search techniques are de- 
scribed in lDib et al.l 1)20091 ). Our best-fit timing solution 
gives a spin frequency v — 0.48277818(5) Hz, frequency 
derivative v = —6.6(1) x 10~ 12 s~ 2 , and second frequency 
derivative v = -6.2(1) x 10~ 18 s" 3 covering the 2008 
Chandra epochs (reference epoch of MJD 54743.0), and 
v = 0.4825962(4) Hz, v = -5.17(5) x 10" 12 s" 2 with neg- 
ligible v for the 2009 Chandra epochs (reference epoch of 
MJD 54854 .0). These values are consistent with those 
reported bv llsrael etail (|2010D and lKaneko et all (|2010f ) 
based on the 2008 Swift and 2009 Fermi observations, 
respectively. 

2.1.2. RXTE Pulsed Fluxes 

Once the timing solution was obtained, we re- 
generated the pulse profiles in the 2-10 keV band, and 
ca lculated the rms pulsed flux according to the formula 
in IDib et alJ (120081 ). using seven harmonics. The results 
are plotted in Figure [T] and reveal a complicated flux 
evolution. In this paper, we focus mainly on time pe- 
riods near the Chandra epochs, and a full analysis of 
the RXTE data will be presented by R. Dib et al. (in 
preparation). Approximately 11 days after the initial 
2008 October 3 (MJD 54742) trigger, during which the 
pulsed flux decayed roughly monotonically, the pulsed 
emission abruptly increased again by 80% between two 
RXTE observations taken on MJDs 54751.2 and 54752.1, 
then decayed monotonically again until it reached a min- 
imum level on MJD 54786.2. As the onset of this second 
event is not resolved, we report an upper limit of 1 day 
for the rise time. This second flux enhancement decayed 
monotonically for ~34 days; the pulsed count rate in 
this period can be parameterized by an exponential fall- 
off with 1/e decay time of ~ 25 days, although a linear 

4 http : //www . atnf . csiro . au/research/pulsar /tempo/ 
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Fig. 1. — Top: rms pulsed flux of IE 1547.0—5408 in the 2-10 keV range obtained with RXTE. The red vertical solid lines mark the 
outbursts in 2008 and 2009, and the arrows indicate the Chandra observation epochs. Bottom: total and rms pulsed Chandra count rates 
in the 1-7 keV band, shown by black squares and red triangles in the upper panels, respectively. The statistical uncertainties are negligible. 
The solid line in 2009 illustrates a power-law fit to the flux decay. Note that the first Chandra observation in 2009 was made with the 
HETG, which precludes a direct comparison of the count rates. Therefore, the data point plotted here is estimated from the spectral 
analysis results. The lower panels show the best-fit blackbody temperature from the PL+BB model, as listed in Table[2] The green dashed 
lines indicate the spin-down rate obtained from the phase-coherent RXTE timing solutions. Uncertainties in v measurements are negligible 
(at 1% level). 



decay is also consistent with the data. Interestingly, the 
2009 event, which is more energetic, is far less dramatic 
as seen by RXTE, with a relatively small increase in the 
pulsed flux observed. 

2.2. Chandra Observations and Results 

The 2008 outburst of IE 1547.0-5408 triggered a series 
of Chandra observations through our Target of Oppor- 
tunity (ToO) program. Five total pointings were made 
on days 4, 7, 15, 18 and 23 after the outburst, with ex- 
posures ranging from 12 to 23 ks. Data were taken with 
the ACIS-S detector in continuous clocking (CC) mode, 
which has a time resolution of 2.85 ms. For the 2009 
event, another Chandra ToO program followed the out- 
burst, with four ACIS-S CC mode observations taken on 
days 2, 4, 8 and 16. The first exposure in 2009 was taken 
with the High Energy Transmission Grating (HETG). 
We include only the zeroth-order events in this latter 
data set in our analysis. A summary of observation pa- 
rameters is provided in Table [TJ 

We performed all the Chandra data reduction using 
CIAO 4.2 with CALDB 4.2.0. We first removed the burst 
intervals in the data. Source counts were extracted from 
a 3" diameter aperture and background counts were from 



the whole chip excluding the central 40" region (i.e., a 
total width of 7'7) . We note that although there is a 
nearby so urce XMMU J155053. 7-5 41925 southwest to 
the pulsar (Gclfand & Gacnslcr 2007), the projected sep- 
aration is always larger than 30" in the CC-mode data. 
Therefore, it does not contaminate the source counts. 
The source count rates in the 1-7 keV energy range are 
reported in Table [T] and plotted in Figure [TJ in which 
the backgrounds have been accounted for, although they 
are less than 0.5%. The count rates are well below 
the telemetry limit of Chandra, and pileup is negligi- 
ble due to the short frame time of the CC-mode expo- 
sures. While the source flux changes between epochs, 
we found no short-term variability within any individ- 
ual exposure. Employing a test algorithm suggested by 
iGreeorv fc Loredol ([l99l . we found < 10% probability of 
source variability within individual observations, which 
is not statistically significant. The flux decay in 2009 can 
be modeled by a power law of index a = —0.21 ± 0.01, 
which is plotted in Figure [T] However, such a simple 
relation is not observed in the 2008 flux evolution. In 
particular, the count rates are very similar between the 
second and third exposures, which is likely related to the 
enhanced pulsed flux detected by RXTE. 
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Fig. 2.— Pulse profiles of IE 1547.0-5408 in l-7keV obtained 
from the Chandra observations, using 64 phase bins. The rms PFs 
from Table [l] are indicated. 

2.2.1. Chandra Timing 

We applied a barycenter correction to the Chandra 
data and then folded the photon arrival times accord- 
ing to the RXTE ephemeris. The resulting lightcurves 
are shown in Figure [2] The pulse profiles in 2008 suggest 
some hints of a multi-peak morphology in the first two 
observations, then evolve into a single peak. By contrast, 
the 2009 profiles exhibit a broad peak at first, with a sec- 
ond peak emerging by the third observation. We found 
no obvious energy dependence of the pulse shape across 
the Chandra band (0.5-10 keV). A direct comparison be- 
tween the profiles in 2008 and 2009 indicates a much 
higher pulse modulation in the former. We estimated 
the rms pulsed fractions (PFs) in 1-7 keV, and the re- 
sults are listed in Table [T] We observe a clear trend of 
increasing pulse modulation as the source recovers after 
the outbursts. 

To look for energy dependence of the modulation, we 
estimated the PFs in the soft (1-3 keV) and hard (3- 
7keV) energy bands separately In the 2008 data, the 
latter shows a systematically higher PF, with a difference 
ranging from APF = 0.04-0.09 (i.e., a 20%-30% change), 
which is statistically significant given the measurement 
uncertainty is only ^0.01. However, such an energy de- 
pendence is not observed in 2009, with the PFs in the 
two bands being consistent with each other. 



Chandra Spectroscopy 
spectra of IE 1547.0-5408 were ex- 



2.2.2. 
The Chandra 

tracted using the tool psextract in CIAO, then binned 
such that every bin has an S/N of at least 10. We 
performed the spectral fits in the 0.5-7 keV range with 
XSPEC vl2.6.0. All nine datasets were fitted jointly with 
a single absorption column density (Ns), We also tried 
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Fig. 3. — Best-fit blackbody plus power-law model fit to the 2008 
and 2009 Chandra spectra (upper panels) with the corresponding 
residuals (lower panels). Different observations are shown by dif- 
ferent colors. The corresponding best-fit spectral parameters are 
listed in Table [2] 

fitting different iV H values for the 2008 and 2009 data, 
and confirmed that they are consistent. We started with 
simple models including an absorbed blackbody (BB) 
and an absorbed power-law (PL), but obtained very poor 
fits (reduced x 2 values over 1.5). An absorbed black- 
body plus power-law (BB+PL) model gives much better 
fits and the results are listed in Table [2] and plotted in 
Figure |3] As shown in the figure, the fit residuals sug- 
gest a hint of a spectral feature ~ 1.3 keV, which is more 
obvious in 2009 than in 2008. However, the significance 
is only ~1<t and deeper exposures are needed to confirm 
this. 

In addition to the BB+PL model, we also considered 
more physical models that account for the Compton up- 
scattering of the thermal photons in the magnetosphere. 
We tried fitting the Resonant Cyclotron Sc attering (RCS; 
iLvutikov fc Gavriill [20061 : iRea et al.ll200ll and the Sur- 
face Thermal Emission and Magnetospheric Scattering 
(STEMS; [Qiell200l IGiiver et al.l [2006T) models to the 
data. In the latter, the gravitational redshift is fixed at 
z = 0.306 during the fit, corresponding to the canon- 
ical neutron star mass of 1.4 Mq and radius of 10 km. 
While these models fit the 2008 data reasonably well, 
the 2009 spectra clearly require an additional hard com- 
ponent. Therefore, we added a PL to the spectral model 
in 2009, with T fixed at 1.33 according to the Suzaku 
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TABLE 2 

Phase-averaged Spectral Parameters of IE 1547.0- 



5408 for Different Models 



Obs. 


7V H 
(10 22 cm" 2 ) 


B 

(10 14 G) 


kT 
(keV) 


r r 

(km) 


/■abs a 77 / jp b 
J ^pl/^th 

(lO^ 11 erg 

s _1 cm -2 ) 










Blackbody - 


- power-law model 







2008 

1 

2 

3 

4 

5 

2009 

6 

7 



4.1±0.1 C 
4.1 ± 0.1 c 
4.1 ± 0.1 c 
4.1 ± 0.1 c 
4.1 ± 0.1 c 

4.1 ± 0.1 c 
4.1 ± 0.1 c 
4.1 ± 0.1 c 
4.1 ± 0.1 c 



0.58 ± 0.03 

0.57 ± 0.02 

0.60 ± 0.04 

55+ 03 
u - dd -0.02 

0.57 ± 0.02 



2.0 ±0.3 

9 1+0.3 

i 7+0.4 
■■••'-0.3 
o 9+0.4 
z - z -0.3 

1.8 ±0.2 



,+0.3 
-0.4 
j+0.3 
-0.4 



2.3] 

2.3^ 
2.8 ±0.3 
2.4 H 
2. 



„+0.3 
-0.5 
=+0.2 
3 -0.3 



0.49 ± 0.03 
0.46 ± 0.02 
0.47 ± 0.02 
44+ - 02 
RCS + hard power-law mode 



4.1 
4.7 
1.0 
4.1 



+0.8 
0.7 
+0.7 
0.5 
+0.6 
0.5 
+0.8 
0.6 



+0.3 
0.4 



2.0] 

2.1 ± 0.2 

2.3 ±0.2 

2.4 ±0.2 



1.92 ±0.03 
1.64 ±0.02 
1.57 ±0.03 
1.55 ±0.03 
1.34 ±0.02 

5.11 ±0.09 
4.61 ±0.05 
3.67 ±0.04 
3.22 ±0.04 



(+1.0 

-0.6 

1+0.8 
-0.5 
,+1.6 
-1.1 

J - z -0.6 
-1.2 
-0.8 



1.4" 1 
l.lj 
2.6 4 



2.1] 



1. 
1. 
2.2 
2.7 



+ 1.8 
-0.8 
+ 1.0 
-0.6 
1.3 
0.8 
+ 1.8 
1.0 



1.08/1342 c 
1.08/1342 c 
1.08/1342 c 
1.08/1342 c 
1.08/1342 c 

1.08/1342 c 
1.08/1342 c 
1.08/1342° 
1.08/1342 c 



2008 

1 

2 

3 

4 

5 

2009 



4.0±0.1 C 
4.0±0.1 C 
4.0±0.1 C 
4.0±0.1 C 
4.0±0.1 C 

4.0±0.1 C 
4.0±0.1 C 
4.0±0.1 C 
4.0±0.1 C 



7+0.05 
-0.04 

0.47 ± 0.04 

1+0.07 
-0.05 



0.47"* 



0.38" 1 



43+ - 07 

u -^°-0.03 
0.41 ± 0.04 



0.32" 1 



,+0.05 
-0.13 

0.17+0-, 16 



-0.04 
,+0.04 
-0.16 
j+0.18 

-0.01 

STEMS + hard power-law model 



0.32] 
0.13"' 



1.33 d 
1.33 d 
1.33 d 
1.33 d 



1.91 ±0.03 
1.63 ±0.02 
1.56 ±0.03 
1.53 ±0.03 
1.33 ±0.02 

5.11 ±0.09 
4.62 ±0.05 
3.67 ±0.03 
3.23 ± 0.04 



0.37 ±0.08 

37+ ' 14 
u -°'-0.12 

30+°- os 

27+ ' 11 
u,z '-0.08 



1.09/1338 c 
1.09/1338 c 
1.09/1338 c 
1.09/1338 c 
1.09/1338 c 

1.09/1338° 
1.09/1338 c 
1.09/1338 c 
1.09/1338 c 



2008 



1 


4.5 ± 0.1 c 


2 61+ - 10 
z '° -0.08 


313 +o.oii 

u.oio_ 00g 


10 d 




1.90 ± 0.03 




1.08/1329 c 


2 


4.5 ± 0.1 c 


2 69+ - 10 
z ' uy -0.07 


n 314+0-OH 

U.D±<i_ oog 


10 d 




1.62 ±0.02 




1.08/1330 c 


3 


4.5 ± 0.1 c 


2.57 ±0.06 


n 90 + O.O2 
U-2»_ g 


10 d 




1.56 ±0.03 




1.08/1331 c 


4 


4.5 ± 0.1 c 


o -+0.2 
°'^-0.3 


Q o.r + 0.05 
U - JO -0.02 


10 d 




1.53 ±0.03 




1.08/1332 c 


5 


4.5 ± 0.1 c 


9 R9 +0.06 
z - DZ -0.04 


n on + 0.005 

u.ou_ 07 


10 d 




1.33 ±0.02 




1.08/1333 c 


2009 


















6 


4.5 ± 0.1 c 


5 4+0-9 
°-*-1.4 


46+ 09 
U - 4D -0.13 


10+ 


1.33+ 


5.08 ± 0.09 


21 + ' 09 
u - zi -0.11 


1.08/1334 c 


7 


4.5 ± 0.1 c 


3 01+ - 12 


0.306 ± 0.003 


10+ 


1.33+ 


4.60 ±0.05 


0.16 ±0.02 


1.08/1335 c 


8 


4.5 ± 0.1 c 


3 6+0-2 


n 019+0.012 
u, ° -0 .002 


10+ 


1.33+ 


3.66 ±0.04 


12+ 003 

u - lz -0.02 


1.08/1336 c 


9 


4.5 ± 0.1 c 


O Q + 0.6 

JU -0.1 


3n+ 02 
"■■'"-con 


10+ 


1.33+ 


3.23 ±0.04 


15+ 002 

u - io -o.oi 


1.08/1337 c 


Note. 


— All uncertainties quoted arc 


90% confidence intervals (i.e. 


I.60-). 









a Absorbed flux in 0.5-7 keV range. 

+* Unabsorbed flux ratio between the power-law and the primary (i.e., blackbody, STEMS, or RCS) components in 0.5-7keV 



range. 

c All nine data sets were fitted jointly with a single column density. 
d Held fixed in the fit. 



results (jEnoto et al.l I2010T ) Fl Compared to the BB+PL 
fit above, these models provide a similar goodness-of-fit 
in terms of the reduced % 2 values. Table [2] lists the key 
parameters of the best-fit models. The scattering optical 
depth t is around 1-2 in 2008 and > 3 in 2009, and the 
thermal velocity j3 of t he electrons is ^0.4-0.5 in 2008 
and - 0.2 in 2009 (see iLvutikov fc GavriiH 120061 for a 
detailed definition of these parameters). 

3. DISCUSSION 

In this paper, we have reported on Chandra observa- 
tions of IE 1547.0-5408 immediately following its 2008 
and 2009 outbursts, along with RXTE timing and pulsed 

5 The INTEGRAL results also suggest that the hard-band PL 
spectral in dex remained stable ove r the period of the Chandra ob- 
servations (den Hartoe ct al. 2009). 



flux behavior following the 2008 outburst and through- 
out the 2009 event. Next we discuss these observations 
in the context of the magnetar model. 

3.1. Spectral and Spin Evolution 

In the twisted magnetosphere model of magnetars 
(|Thompson et al.1 [2002) , the observed X-ray luminos- 
ity of a magnetar is determined both by its surface 
temperature and by magnetospheric currents, the lat- 
ter due to the twisted dipolar field structure. The sur- 
face temperature in turn is determined by the energy 
output from within the star due to magnetic field de- 
cay, as well as on the nature of the atmosphere and 
the stellar magnetic field strength. This surface ther- 
mal emission is resonantly scattered by the current par- 
ticles, thus resulting in an overall spectrum similar to a 
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Comptoni zed bl ackbodv (e.g. ILvutikov fe Gavriiil 120061 : 
iRea et al.l 120081 iZane et all |2009T>. In this model, the 
greater the twist angle, the greater the scattering, the 
harder the spectrum, and the greater the X-ray lumi- 
nosity Lx- In addition, the surface heating by return 
currents is believed to contribute substantially to Lx, 
at least at the same level as the the rmal component 
induc ed from the interior field decay (| Thompson et al] 
2002). Magnetar outbursts in this picture occur with 
sudden increases in twist angle, consistent with the 
gene ric hardening of m a gnetar spectra during outbursts 
(e^ lKaspi eTaIT[200l 1 Woods et all [200l llsrael et al] 
[20071) . 

Other observational evidence provided in support 
of the twisted magnet osphere model as proposed by 
iThompson etHI ([2002) is a correlation between mag- 
netar spectral hardness and spin-down-inferred mag- 
netic field strength B, when comparing different so urces 
(IMarsden fe White! 12001 IKaspi fc BovdsUml [2010h . In 
this case, B is an observational proxy for the magne- 
tospheric twist angle. On the other hand, some mag- 
netars have shown dramatic spin-down rate variations, 
with order-of-magni tude changes in v seen on a vari- 
ety o f timescales (e.g. lGavriil fc K aspi 2004; Wood s~et al] 
2004). The origin of these variations is unknown. Nev- 
ertheless, in the context of the twisted magnetosphere 
model, a varying twist angle might naively be expected 
to be accompanied by a changing v (due to chang- 
ing effective B), and corresponding spectral and flux 
changes. However, some decoupling between v and the 
radiative behavior might be expected, particularly as the 
spin-down is affected most by a narrow field-line bun- 
dle near the light cylinder, whereas the radiation orig- 
inates from the surface. Field-line tw ists likely prop- 
agate outward (jThompson et al.ll2002fl . which suggests 
the radiative changes should oc cur prior to v changes 
(jBeloborodov fc Thompson1[2007l ). 

In contrast to the picture in which a magnetar outburst 
is accompanied by an enhanced magnetospheric twist, 
IQzel fc Giiverl ([2007D suggest that in outburst, the mag- 
netosphere may be stable, with radiative evolution being 
due to changes in the surface thermal emission. Using 
a spectral model consisting of a resonant Comptonized 
atmosphere- modified blackbody (the STEMS mod el), 
fits to data for XTE J1810-197 (|Guver et alj|2007l ) re- 
sult in the spectrally inferred B being stable, with all 
radiative changes being due to changes in the surface 
thermal emission. 

For IE 1547.0—5408, the source spectrum only showed 
significant variability over a short perio d of time (~1 day) 
following the 2008 and 2009 outbursts ([Israel et alj|20l'ot 
Scholz et al., in preparation), but remained stable over 
the Chandra observations, during which the flux changed 
substantially (see Table [2]) . The latter seems opposite to 
the hardness/flux correlation predicted by the twisted 
magnetosphere model. On the other hand, we note that 
the flux decay over the Chandra epochs has a compara- 
ble timescale to that of the pulse profile variations (Fig- 
ure [5]) , which generally agrees with the predictions of the 
magnetar model. More intriguing are the timing results. 
As described in Section [2. 1.11 we found a factor of 2.2 in- 
crease in the spin-down rate \ v\ between the first and last 
Chandra observations in 2008, a substantial change even 



by magnetar standards. For purely dipole spin-down, 
this naively implies a ~50% increase in the effective sur- 
face dipole field strength, from an initial spin-inferred 
value of 2.8 x 10 14 G at the epoch of the first Chandra 
observation, to a value of 4.1 x 10 14 G at the epoch of 
the last. This is contrary to SGR 1806—20, in which 
the spectral response lagged behind the tor que variation, 
sugge sting some hysteresis in the system (|Woods et al I 
2007) . In our case, the lack of associated spectral changes 
is unexpected in the twisted magnetosphere model, un- 
less the spin-down is decoupled from the site of t he ra - 
diative events, as sugg ested b y IGavriil fc Kaspil (12004) 
for IE 1048.1-5937 (see also iBeloborodovl I2009D . Al- 
though B fields inferred from long-term spin-down have 
been used to comp are with those inferred spectrally (e.g. 
IGiiver et al.H2007f ). our results call into question the reli- 
ability of such a comparison when using short-term spin- 
down rate. 

Indeed, for the 2008 Chandra observations of 
IE 1547.0-5408, spectral fits using the STEMS model 
yield a very similar B field value for all the observa- 
tions (see Table ^ in spite of the strongly varying v. 
One way to reconcile this is to interpret the B field 
measured spectroscopically as higher-order multipoles in 
localized X-ray-emitting regions rather than the global 
dipole field responsible for spin-down. However, we note 
that the B field obtained from the STEMS model is 
lower than the spin-down-inferred value at the epoch of 
the last 2008 Chandra observation, which is difficult to 
be explained by the above picture. As an alternative, 
the extra s pin-down torque coul d be attributed to parti- 
cle winds ( Harding^eFal] [1999) . Compar ing the spin- 
down rates of IE 1547.0-540 8 in 2007 (jCamilo et al] 
I2008D . 2008 and 2009 (Section |2~TT|) . it is obvious that 
the spin-down torque changed drastically between these 
epochs. Torque variations have been o bserved in magne- 
tars and ordinar y rad i o pulsars (e.g. IKaspi et al.l [20031: 
Gavriil fc Kaspil [20041: Kramer et al] 120061: iLvne et al] 



20101 : iLivingstone et al. These could be related 



to changes in plasma conditions in the magnetosphere, 
which may not necessarily have any observable effects in 
th e radiative pro p erties (see ?). Based on Equation (9) 
in Hardin g" et al.l |l999), a factor of 2.2 increase in \i>\, 
and hence in E, requires a steady wind luminosity (L p ) 
1.5 times larger than the dipole spin-down luminosity, 
implying L p k 1.5 x 10 35 ergss -1 . For a typical X-ray 
efficiency of below 1%, the particle-induced flux would be 
less than 10~ 12 ergscm^ 2 s _1 , much smaller when com- 
pared to the source's flux (see Tabled]), thus, unlikely to 
be detected. 

In any case, the absence of spectral variations in the 
presence of flux changes remains puzzling whether or not 
the magnetospheric twist angle varied in the outburst. 
Moreover, a reasonable model will need to explain why i> 
changed drastically following the 2008 event, but stayed 
constant in 2009. This may reinforce the requirement of 
decoupling between the spin-down and the source of the 
radiative changes, hence presumably between the spin- 
down-inferred magnetic field and that inferred spectrally. 

3.2. Flux Evolution During the 2008 and 2009 Events 

The 2008 and 2009 events exhibited very different flux 
evolutions. Immediately after the 2009 outburst, the per- 
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Fig. 4.— rms PF vs. absorbed flux of IE 1547.0-5408 in the 
1-7 keV band. The filled triangles and squares show the 2008 and 
2009 Chandra observations, respectively. The green open circle and 
red do ts represent results f rom XMM-Newton and Swift, respec- 
tively (Halpcrn ct al. 2008; S. A. Olausen et al., in preparation; 
P. Scholz et al., in preparation), with the XMM-Newton epochs 
labeled. Note that the 2007 and 2008 XMM-Newton data points 
are only upper limits, since the values reported in the literature are 
area PF estimates, which are larger than the rms PFs by a con- 
stant which depends on the pulse shape. The dashed line shows a 
power-law fit to the Chandra data. 

sistent flux increased by a factor of ~ 500 (Scholz et al. 
in preparation), while the pulsed flux evidently showed 
only very little variation (less than a factor of two). We 
found a monotonic flux decay during the 2009 recovery, 
with a power-law of index —0.21 ± 0.01, which is corn- 
arable to -0 306 ± 0.005 for CXOU J164710.2-455216 
Woods et aL l lgOlin. but not a s steep as -0.69 ± 0.03 for 
IE 2259+58 6 dZhu et al.l[2008h or -0.92 ± 0.02 for XTE 
J1810-197 HWoods et al.ll2005ft . 

In contrast, the 2008 event is less energ etic; the to- 
tal fl ux increased only by a factor of ~100 ([Israel et al.l 
I20l0t) . but our RXTE results reveal a pulsed flux varia- 
tion by a factor of ~4, far greater than in the 2009 event. 
Also, the flux decay in 2008 showed a more complicated 
history. As is clear in Figure [1] we observed an addi- 
tional flux enhancement around MJD 54752, ~11 days 
a fter the initial trigger, R l asting for ~30 days. 

lEichler fc Sha isultanov (2010) suggest that radiative 
outbursts in magnetars could generally be preceded by 
glitches, with the delay between the two events depend- 
ing on the depth at which the glitch-induced energy re- 
lease occurs. For the initial event in 2008, we are un- 
able to tell whether a glitch preceded the radiative out- 
burst, due to the lack of RXTE observations prior to the 
outburst. However, we can rule out any glitch between 
the initial 2008 event and the second flux enhancement 
10 days later. It is possible that the initial event actu- 
ally involved glitches occurring in two different places in 
the stellar interior, at substantially different depths, such 

6 Although th is second enhancement is not reported by 
llsrael et al.1 1120101) in their study of the Swift observations cov- 
ering the same period, we have re-analysed the same Swift datasct 
and confirmed the pulsed flux enhancement we observe with RXTE 
and Chandra (P. Scholz et al., in preparation). 



that the delays from the glitch to the X-ray enhancement 
were different, but we note that this picture does not ex- 
plain the sharp rise of the second flux enhancement. 

3.3. Pulsed Fraction Evolution 

Our results in Section [2] clearly indicate a strong anti- 
correlation between the PF and the phase-averaged X- 
ray flux, at least during the 2008 and 2009 outbursts. 
This is plotted in Figure IH and suggests an approxi- 
mate power-law relation between the two observables. 
The trend is also supported by the XMM-Newton and 
Swift measurements taken in the same period (S. A. 
Olausen et al., in preparation, P. Scholz et al., in 
preparation). Similar anti-correlations have been ob- 
served in IE 1048.1-5937 and CXO U J164710.2-455216 
(|Tam et alj|2008t llsrael et al.ll2007[ ). while positive corre- 
lations were found in XTE J1810-197 and IE 2259+586 
(|Gotthelf et al.l l200l iZhu et al.l 12051) . This variety of 
behaviors is consistent with the picture in which, from 
source to source, the location and geometry of the re- 
gion on the star affected in the outburst are different. 
Previous studies proposed that an anti-correlation could 
be the consequence of an increased emitting area due to 
an outburst, such that part of the hot spot becomes vis- 
ible at any phase, thus r educing the pulse modulation 
(e.g. lHalpern et all 120081 ). We note that this scenario 
depends critically on the location of emission zone on 
the stellar surface as well as the viewing geometry (see 
iBogdanov et al.ll2008f) ; it may be possible to obtain ei- 
ther a monotonic increase or decrease of PF for the same 
area 'hot spot' depending on its location on the stellar 
surface. 

We point out that the trend of decreasing PF for in- 
creasing phase-averaged X-ray flux observed in the Chan- 
dra data does not seem to hold for other observations of 
IE 1547.0 -5408 before 2 08. B ased on XMM-Newton 
exposures, lHalpern et al.l (|2008) reported area PFs of 
15% in quiescence in 2006 and 7% in the high state in 
2007. As shown in Figure HJ these values deviate sig- 
nificantly from the 2008-2009 trend. The discrepancy 
seems too large to be reconciled by a difference in the in- 
strument response, and is even larger if the differing PF 
estimate methods are accounted for. It is possible that 
the outburst in 2008 induced some permanent changes 
in the i3-field configur ation or emission geo metry. The 
pulse profiles shown in lHalpern et al.l (|2008[ ) also appear 
to have a different shape from the ones shown in Figure^ 
providing further support to this picture. 

4. CONCLUSIONS 

We have presented results from Chandra and RXTE 
observations of IE 1547.0-5408 following its 2008 and 
2009 outbursts. These allow a direct comparison between 
the two events. We found that over the 2008 Chandra 
observation epochs, the pulsar spin-down rate increased 
by a factor 2.2, in the absence of corresponding spec- 
tral changes, whereas such variation in v is not observed 
after the more energetic 2009 event. This provides evi- 
dence of decoupling between magnetar spin and radiative 
properties. The absence of spectral changes simultane- 
ous with significant flux decay is surprising for models of 
magnetar outbursts. Our results also revealed a strong 
anti-correlation between the PF and phase-averaged flux 
of the source. While both 2008 and 2009 data follow 



Ng et al. 



the same trend, pre-2008 measurements show significant 
deviation, suggesting that the 2008 outburst may have 
induced permanent changes in the emission geometry. 
Finally, we note that IE 1547.0—5408 demonstrated sig- 
nificant spectral changes only within the first day after 
the 2008 and 2009 events, which highlights the impor- 
tance of prompt observations in future studies for under- 
standing post-outburst relaxation of magnetars. 
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